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Abstract. 

We extend the standard fiber bundle model (FBM) with the local load sharing in such a 
way that the conservation of the total load is relaxed when an isolated fiber is broken. In 
this modified FBM in one dimension (ID), it is revealed that the model exhibits a well-defined 
phase transition at a finite nonzero value of the load, which is in contrast to the standard ID 
FBM. The modified FBM defined in the Watts- Strogatz network is also investigated, and found 
is the existences of two distinct transitions: one discontinuous and the other continuous. The 
effects of the long-range shortcuts are also discussed. 



The fiber bundle model (FBM) was introduced to explain a variety of phenomena related 
with the problems of damage spread [1]. More specifically, consider a system composed of a 
bundle of mingled fibers exposed to a common stress along the direction of fibers. It is natural 
to assume that each fiber has a certain value of the capacity, or the threshold value of the load, 
and if the load on the fiber exceeds the capacity allowed to it, the fiber is broken. Since the 
broken fiber stops functioning, its neighbor fibers should carry the load which was assigned to 
the broken one [2] . The damage caused by the breaking of a single fiber may propagate over 
its local neighbors and then stop spreading. Another interesting possibility is the existence 
of an avalanche or a cascading failure, i.e., the original local damage spreads over the whole 
system, causing the complete breakdown of the total fibers. Recent blackout in the United 
States can also be understood as the cascading failures of power cables and power stations 
caused by the spread of the overload breakdowns of elements in power grid. 

Recently, the subject of complex networks is the one of the most studied research area, 
not only for its intrinsic interest but also for its promising broad future applicability [3]. 
When a dynamic system is defined on a complex network, the emergent behavior may totally 
differ from the one on the local regular array which has been the standard arena of statistical 
physics. For example, the phase transitions of the Ising and the XY models on complex 
networks have been shown to belong to the mean-field universality class [4]. The spread of 
damage on networks, which is not described by the equilibrium Hamiltonian, also draws much 
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interest in relation to the epidemic spread [5] as well as the overload breakdowns of networks 
subject to local [6] or nonlocal [7] loads. 

The fiber bundle model in this work on a general network structure is defined as follows: 

1. To each vertex in the network of the size N the capacity c, for the ith vertex is assigned 
at random following a given distribution function (in most cases we use the uniform 
distribution with c$ € [0, 1]). 

2. The system is then exposed to a total load Na, and thus each vertex should carry the 
same load Oj = a (we in this work a — is used as an initial condition). 

3. Examine the breaking condition <n > c; for each i and if the inequality is satisfied the 
vertex i is broken and its load <Ji is equally transfered to its surviving neighbor vertices. 
For example, if the vertex i has ki surviving neighbors the vertex j, which is a neighbor 
of z, has now new load: <7j —> aj + ai/ki. Repeat this step until there is no more 
breaking. 

4. Increase the total load Na — > N(a + da). The total load should be carried by unbroken 
surviving vertices, the number of which is N s , and accordingly each vertex is assigned 
the new load ai — * Oj + (N/N s )da. Note that we denote the local load on the vertex i 
as (Ji while a is the total load divided by N. 

The above procedures are identical to the standard FBM with the local load sharing (LLS) 
in the literature [2]. However, our modified FBM in the present work differs as follows: 

• When an isolated vertex, which has no surviving neighbor vertices, is overload broken, 
we assume that the load which has been carried by this vertex is not transferred to 
other vertices. On the other hand, in the standard FBM with LLS the load transfers 
via edges, not via chain of surviving vertices, and consequently the load of the broken 
isolated fiber can still be transferred to other (suitably defined) closely located vertices. 

Our modified FBM appears to be reasonable in itself since the standard FBM assumes that 
the load can be transferred to vertices which are not connected to the broken vertex through 
a chain of surviving vertices. In contrast, our model allows only the load transfer across the 
actually existing surviving vertices. One disadvantage of our modified FBM is that the total 
load applied to the system is not constant in time. It might be difficult to find the actual 
physical realization of our model in the original context of the FBM. However, in abstract 
networks such as Internet one may not rule out the possibility that when an isolated vertex 
is broken its load does not spread to other vertices but just disappears. It should be noted 
that even in our modified FBM, the total load is resumed when each vertex is assigned the 
new load as <Ji — > <7, + (N/N s )da (see step 4) since the breaking of an isolated vertex also 
decreases N s . Accordingly, the total load conservation condition is relaxed only during the 
spread of damage in step 3. 

We repeat the above described procedure for our modified version of FBM and measure 
the ratio p = N s /N as a function of the load a averaged over 10 4 ~ 10 5 different capacity 
distributions and different network structures. In most numerical simulations, we use da = 
0.001, which has been numerically confirmed to be small enough, and network sizes are N = 
128, 256, • • • , 524 288 (or from N = 2 7 to 2 19 ). 

We first consider the local regular one-dimensional (ID) network with the connection 
range r. For various values of r, including r = 1 which corresponds to the usual ID lattice, 
we observe the qualitatively the same behavior. In Fig. for r = 2, it is shown that 
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Fig. 1 Fig. 2 

Fig. 1 - The modified fiber bundle model on the ID regular lattice with couplings up to the next 
nearest neighbors with the uniform distribution of the capacity a 6 [0, 1]. (a) Density p of unbro- 
ken fibers versus the load a for N = 2 7 , 2 8 , • • • , 2 19 (from top to bottom). As a is increased, the 
network begins to break up, resulting in p = at sufficiently large values of a. (b) p versus N for 
a = 0.2812,0.2814,0.2816,0.2818, and 0.2820 (from top to bottom). Around a w 0.2817, p decays 
algebraically with N, manifesting the existence of the continuous phase transition, (c) Precise deter- 
mination of Oc from the finite-size scaling form Determined are a = 0.50(1) and er c = 0.28165(2). 
(Curves correspond to N = 2 12 , 2 13 , • • • , 2 19 , respectively, from bottom to top on the left-hand side 
of the crossing point), (d) Scaling collapse of all data points in (c) through the use of the finite-size 
scaling v = 2.0(1) is obtained. 

Fig. 2 - The FBM on the ID regular lattice with connections up to the next nearest neighbors with 
the Gaussian distribution of the capacity (the mean and the variance are fixed the same as in Fig. ^ . 
Finite-size scaling (compare with Fig. for the uniform distribution) yields a c — 0.27363(2) and 
critical exponents a = [3/v = 0.50(1) and v — 2.0(1), identical to the uniform distribution. 



at a certain value of the load a the ratio p between the number of working fibers and the 
total number of fibers vanishes beyond some critical value a c . As the system size N becomes 
larger [N = 2 7 , 2 8 , • • ■ , 2 19 (— 524 288) in Fig. da)], the phase transition becomes clearer. It 
is interesting to note that the existence of the phase transition at nonzero a c is completely 
different from the well-established result of a c — for the standard FBM in the thermodynamic 
limit [2] . Simple change of the conservation law of the total load yields the dramatic change 
of the universality class. 

We in Fig. ^ (b) show p as a function of N in log scales at various values of a = 
0.2812,0.2814,0.2816,0.2818, and 0.2820 (from top to bottom). For a > a c w 0.2817, p(N) 
shows a downward curvature, implying p — > in the thermodynamic limit, while p(N) bends 
upward for a < a c , suggesting that p(N) saturates to a finite value as N becomes larger. 
Separating these behaviors, there exists a well-defined a c at which p as a function of N decays 
algebraically, and the system does not possess any length scale except N at a = a c . This is 
commonly found behavior in the system with a continuous phase transition. 

For a more precise determination of er c , we use the standard finite-size scaling method with 
the scaling assumption: 

p(a,N)=N- a f({a-a c )N 1 ^), (1) 

where f(x) is the scaling function and the exponent v describes the divergence of the cor- 
relation length £ [or more generally, the correlation volume for the case when the length is 
ill-defined (see discussions in, e.g., Ref. [4])] near the critical point, i.e., £ ~ \a — <r c \~ v . From 
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Fig. 3 - The FBM on the Watts-Strogatz small-world networks with the connection range r — 2 
and at various rewiring probabilities P — (a) 0.04, (b) 0.1, (c) 0.2, and (d) 0.5. Except (a), all 
show interesting transition behaviors: the density of working fibers p drops down abruptly at a 
certain value of the load <7i and then approaches p = continuously at 02. Networks with the sizes 
TV = 2 7 , 2 8 , • • • 2 19 (from top to bottom in each plot) and the uniform capacity distribution are used. 



the definition of the critical exponent /?, p ~ (<r c — a) 13 for a < a c in the thermodynamic limit, 
we obtain the relation a — (3 jv. Figure ^c) is for a precise determination of a c from the 
finite-size scaling (QJ: pN a is shown to have a unique crossing point at a = a c = 0.28165(2) 
with the exponent a = (3/v — 0.50(1). Once <r c and the exponent a are determined, one can 
again use the finite-size scaling £Q) to make the all data points in Fig. UH C ) collapse to a single 
smooth curve as displayed in Fig. ^d) with the determination of v = 2.0(1). 

The above results of the existence of the continuous phase transition with critical exponents 
v k 2 and j3 « 1 is expected to be universal regardless of the detailed form of the distribution 
function of the capacity. We in Fig. |2 show the finite-size scaling plot for the same FBM but 
with the Gaussian distribution of the capacity instead. For comparisons, the same values of the 
mean and the variance as ones with the uniform distribution in Fig. ^ are used. Confirmed is 
the existence of the continuous phase transition with the identical critical exponents. Although 
we have not used the commonly used Weibull distribution, we believe that the universality 
class should be identical for the various capacity distribution functions. 

We next study the FBM on the Watts-Strogatz (WS) networks, which is generated fol- 
lowing the standard procedure in Ref. [8]. The two important parameters in the WS network 
are the connection range r, (r = 2 is used here) and the rewiring probability P controlling 
the number of long-range shortcuts. As soon as P has a nonzero value, it is well-known that 
the network undergoes so-called small-world transition that the network diameter increases 
logarithmically with the network size N. In Fig.^p versus a is shown at P = (a) 0.04, (b) 
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Fig. 4 - The FBM on the WS network with the rewiring probability P = 0.5. (a) The maximum of 
the derivative of p(a) with respect to a versus the system size N. The divergence of the derivative 
implies the discontinuous transition nature in the thermodynamic limit. a\ is defined as the value of 
a when this maximum occurs, (b) The location o\ of the derivative maximum versus 1/N. In the 
thermodynamic limit d ~ 0.207 is estimated. 

Fig. 5 - The continuous phase transition in the FBM on the WS network with P — 0.5. (a) 
pN a versus a at various sizes N = 2 13 ,---2 19 with a = 0.5 crosses at a unique crossing point, 
resulting in the estimation <T2 = 0.22376(2) as the critical point of the second phase transition of the 
continuous nature, (b) Finite-size scaling collapse Q of data points in (a). The critical exponents 
a = J3/v = 0.5, v = 2.0 (and thus f3 = 1.0) result in the scaling collapse. 



0.1, (c) 0.2, and (d) 0.5. At sufficiently large values of P, it is clearly demonstrated that the 
system has two phase transitions at o\ and 02 (01 < At a\, p appears to drop abruptly 
to a nonzero finite value and as a is further increased, p eventually vanishes beyond oi- 

To investigate the nature of the first transition at o\, we compute dp/ da as a function of 
a and then measure at which value of a the derivative takes the maximum value. In Fig. Ufa), 
the maximum value of the derivative (dp/da)^^ is plotted as a function of the system size 
N: The power-law divergence (dp/da) max ~ N 0A clearly suggests that the transition nature 
is indeed discontinuous in the thermodynamic limit of N — * 00. The location of the derivative 
maximum, a±, is displayed as a function of 1/N in Fig. Ufb); ci ~ 0.207 is concluded in the 
thermodynamic limit. We then apply the finite-size scaling analysis to the regime where p 
approaches zero continuously (see Fig.EJ. The unique crossing point in the pN a versus a at 
various N in Fig.^Ia) as well as the smooth scaling collapse in Fig.OJb) is again observed with 
the identical values of critical exponents, i.e, a w 0.5, v w 2.0, and (3 « 1.0. It is noteworthy 
that the FBM on the scale-free network of Barabasi and Albert (BA) [9] has shown to have 
a different transition nature [6], i.e., a discontinuous phase transition without the subsequent 
continuous transition as in the WS network considered here. Consequently, as soon as the 
network undergoes the global spread of the damage, no fiber survives the disaster in the 
BA network, while in the WS network fibers which survived the discontinuous transition are 
subject to the second continuous transition at a higher load. 

In the context of complex networks, the size of giant component has been more frequently 
measured since it can detect the global breakdown of the information flow between two ar- 
bitrarily chosen vertices. In the viewpoint of the universality, the use of the different order 
parameters is expected to change neither the nature of the phase transition nor the values of 
critical exponents, which are the main interest in this work. The density of surviving fibers 
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in the WS network can also have a practical meaning if one considers the real fiber bundles 
in which most inter-fiber couplings are local while some of them are long-ranged. 

Figure [B] summarizes the phase diagram obtained in the present work. There exist in 
general two phase transition lines separating three different phases (two unbroken phases UI 
and UII with p ^ and the broken phase B with p — 0): <j\ is estimated from the maximum 
of dp I da while oi is obtained from the standard finite-size scaling analysis applied for the 
continuous phase transition. As the load a is increased from zero, the system only suffers 
from locally isolated damage spreading (the phase UI in Fig. |HJ) until the first discontinuous 
transition point <j\ is reached. As a crosses <j\ from below, the system undergoes abrupt 
change, probably in the form of the cascading failures or an avalanche, and the number of 
fibers still alive decays discontinuously. Interestingly, the existence of the long-range shortcuts 
not only makes the discontinuous transition happen, but also help the system to avoid complete 
disaster, i.e, the global spread of damage does not destroy all fibers in the system but there 
still remain a finite fraction of fibers which survived the first transition at a\ (the phase UII in 
Fig. El • Further increase of the load then eventually destroys all fibers resulting in the p = 
state (phase B in Fig. EJ. At sufficiently small values of P, it appears that the system has 
only one phase transition which is continuous in nature: At P < 0.04, it is found that the 
maximum of dp/ da and the continuous phase transition occur simultaneously. At P > 0.06, 
the system starts to have well-separated two phase transitions, and the increase of P appears 
to make the system more strong (both o\ and 02 increase with P in the intermediate region) . 

In summary, we in this work have extended the existing model of fiber bundles in such 
a way that the total load conservation condition is relaxed. For a pure ID regular network, 
our modified FBM has been found to show well-defined continuous phase transition at a finite 
nonzero value of the load a c with the critical exponents (3 ~ 1.0 and v w 2.0. This is in a 
sharp contrast to the standard FBM where a c — has been confirmed in the thermodynamic 
limit. It has also been verified that the universality class of the phase transition does not 
change when different distribution function for the capacity is tried. For the Watts-Strogatz 
network, as the rewiring probability P is increased, the system undergoes two phase transitions 
at a\ and 02 • At o\ the density p of working fibers decreases discontinuously to a nonzero 
finite value. At <Ji(> cti), p approaches continuously zero, with the phase transition with 
the identical universality class (/3 « 1.0 and v ss 2.0). The roles played by the existence of 
long-range shortcuts have been shown to be two-fold: On the one hand, long-range shortcuts 
make the fiber breaking transition abrupt, helping the damage spread across the system more 
easily. On the other hand, the long-range shortcuts prohibit the system from being completely 
destroyed, and even after the global propagation of damage (at <7i), there still remain a finite 
fraction of fibers still working. 

The author thanks D.H. Kim for useful discussions. This work has been supported by 
the Korea Science and Engineering Foundation through Grant No. R14-2002-062-01000-0 and 
Hwang-Pil-Sang research fund in Ajou University. Numerical works have been performed on 
the computer cluster Iceberg at Ajou University. 
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